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Abstract: The high parton density effects are strongly dependent of 
the spatial gluon distribution within the proton, with radius R, which cannot 
be derived from perturbative QCD. In this paper we assume that the unitarity 
corrections are present in the HERA kinematical region and constrain the 
value of R using the data for the proton structure function and its slope. 
We obtain that the gluons are not distributed uniformly in the whole proton 
disc, but behave as concentrated in smaller regions. 
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The description of the dynamics at high density parton regime is one of 
the main open questions of the strong interactions theory. While in the region 
of moderate Bjorken x (x > 10~ 2 ) the well-established methods of operator 
product expansion and renormalization group equations have been applied 
successfully, the small x region still lacks a consistent theoretical framework 
(For a review see Basically, it is questionable the use of the DGLAP 
equations @, which reflects the dynamics at moderate x, in the region of 
small values of x, where the gluon distribution determines the behavior of 
the observables. The traditional procedure of using the DGLAP equations 
to calculate the gluon distribution at small x and large momentum trans- 
fer Q 2 is by summing the leading powers of ot s lnQ 2 In(-), where a s is the 
strong coupling constant, known as the double-leading-logarithm approxima- 
tion (DLLA). In axial gauges, these leading double logarithms are generated 
by ladder diagrams in which the emitted gluons have strongly ordered trans- 
verse momenta, as well as strongly ordered longitudinal momenta. Therefore 
the DGLAP must breakdown at small values of x, firstly because this frame- 
work does not account for the contributions to the cross section which are 
leading in a s ln(-) @. Secondly, because the parton densities become large 
and there is need to develop a high density formulation of QCD H, where 
the unitarity corrections are considered. 

There has been intense debate on to which extent non-conventional QCD 
evolution is required by the deep inelastic ep HERA data |l| . Good fits to the 
F 2 data for Q 2 > 1 GeV 2 can be obtained from distinct approaches, which 
consider DGLAP and/or BFKL evolution equations || H|. In particular, the 
conventional perturbative QCD approach is very successful in describing the 
main features of HERA data and, hence, the signal of a new QCD dynamics 
has been in general hidden or mimicked by a strong background of conven- 
tional QCD evolution. For instance, recently the magnitude of the higher 
twist terms was demonstrated to be large in the transverse and longitudinal 
structure function, but as these contributions have opposite signal the effect 
in the behavior of F 2 structure function is small [0]. At this moment there 
are some possible signs of the high density dynamics in the HERA kinemat- 
ical region: the behavior of the slope of structure function and the energy 
dependence of diffractive structure functions ||. However, more studies and 
precise data are still needed. 

Our goal in this letter is by assuming the presence of the high density 
effects in the HERA data, to constrain the spatial distribution of the gluons 
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inside the proton. The radius R is a phenomenological parameter which is 
not present in the linear dynamics (DGLAP/BFKL) and is introduced when 
the unitarity corrections are estimated. In general, the evolution equations at 
high density QCD (hdQCD) |4| resum the powers of the function n(x, Q 2 ) = 
2qTr5 %G(x, Q 2 ), which represents the probability of gluon-gluon interaction 
inside the parton cascade. At this moment, these evolution equations 

• match the DLA limit of the DGLAP evolution equation in the limit of 

low parton densities (k — > 0); 

• match the GLR equation in first order of unitarity corretions [(9(k 2 )]. 

Although the complete demonstration of the equivalence between these for- 
mulations in the region of large k is still an open question, some steps in 
this direction were given recently ||, W\. One of the main characteristics 



in common of these approaches is the behavior of the structure function in 
the asymptotic regime of very high density [K|: F2(x,Q 2 ) oc Q 2 R 2 ln(l/x). 



Therefore, although the parameter R cannot be derived from perturbative 
QCD, the unitarity corrections crucially depend on its value. Here we will 
discriminate the range of possible values of R considering the AGL approach 
for the high density systems and the HERA data for the structure function 
F 2 (x, Q 2 ) and its slope. In the HERA kinematical region the solution of the 
AGL equation is identical to the GLR solution [[11]], which implies that our 
estimates in principle are not model dependent. 

We will consider initially the physical interpretation of the R parameter, 
and present later a brief review of the AGL approach and our estimates. 
The value of R is associated with the coupling of the gluon ladders with the 
proton, or to put it in another way, on how the gluons are distributed within 
the proton. R may be of the order of the proton radius if the gluons are 
distributed uniformly in the whole proton disc or much smaller if the gluons 
are concentrated, i.e. if the gluons in the proton are confined in a disc with 
smaller radius than the size of the proton . In a first approximation, the 



radius is expected to be smaller than the proton radius. This affirmative is 
easy to understand. Consider the first order contribution to the unitarity 
corrections presented in Fig. [l], where two ladders couple to the proton. The 
ladders may be attached to different constituents of the proton or to the 
same constituent. In the first case [Fig. [I] (a)] the unitarity corrections are 
controlled by the proton radius, while in the second case [Fig. [1] (b)] these 
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corrections are controlled by the constituent radius, which is smaller than 
the proton radius. Therefore, on the average, we expect that the radius will 
be smaller than the proton radius. The value of R 2 reflects the integration 
over b t in the first diagrams for the unitarity corrections. 

In our estimates for the parameter R we will use the AGL approach |L1 



Here we present only a brief review of this approach and refer the original 
papers for details. Basically, in the AGL approach the behavior of the gluon 
distribution can be obtained from the cross section for the interaction of a 
virtual gluon with a proton. In the target rest frame the virtual gluon at 
high energy (small x) decay into a gluon-gluon pair long before the interaction 
with the target. The gg pair subsequently interacts with the target, with the 
transverse distance r t between the gluons assumed fixed. It allows to factorize 
the total cross section between the wave function of the virtual gluon and the 
interaction cross section of the gluon-gluon pair with the target. The gluon 
wave function is calculable and the interaction cross section is modelled. 
Considering the unitarity corrections for the interaction cross section results 



that the gluon distribution is given by the Glauber-Mueller formula |TT 

xG(x, Q 2 ) = A r \ ^- < -^{C + ln(K G (x', - 2 )) + EkfaV, - 2 ))} , (1) 
7T Jx x J irr t r t r t 

where C is the Euler constant, E\ is the exponential function and the function 
i^g{ x i -pi) = fUS" 71 " 7 "* xG(x, -7). If equation ([[]) is expanded for small kq, the 
first term (Born term) will correspond to the usual DGLAP equation in the 
small x region, while the other terms will take into account the unitarity 
corrections. The Glauber-Mueller formula is a particular case of the AGL 



equation proposed in Ref. flip , and a good approximation for the solutions 
of this equation in the HERA kinematical region, which we will use in this 
work. 

A similar procedure can be used to estimate the unitarity corrections for 
the F2 structure function and its slope. In the target rest frame the proton 
structure function is given by |14| 

F 2 (x,Q 2 ) = J?— fdz f d 2 r t \^(z,r t )\ 2 ^(z,r t ) , (2) 

where \& is the photon wave function and the cross section a q ^(z, r 2 ) describes 
the interaction of the pair with the target. Considering only light quarks 
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(z = u, d, s), the expression for \1/ derived in Ref. |TJj and the unitarity 



corrections for the interaction cross section of the qq with the proton, the F2 
structure function can be written in the AGL approach as [TT1 



^(*,Q 2 ) = ^ E 4 j}^T{C + ln(K q (xy t )) + E 1 ( Kq (x^))} t (3) 

'" f=u,d,s cp * 



where the function n q (x,r 2 ) = 4/9kg(%, rf). The slope of F 2 structure func- 
tion in this approach is straightforward from the expression (|3|). We obtain 
that 

' IF2U - Q2] - R2Q2 Y,e){C + ln^ q {x 1 Q 2 )) + E l ^ q {x,Q 2 ))} . (4) 



dlogQ 2 2tt 2 



u,d,s 



The expressions and (|j) predict the behavior of the unitarity correc- 
tions to F2 and its slope considering the AGL approach for the interaction of 
the qq with the target. In this case we are calculating the corrections asso- 
ciated with the crossing of the qq pair through the target. Following |16J we 
will denote this contribution as the quark sector contribution to the unitarity 
corrections. However, the behavior of F2 and its slope are associated with the 
behavior of the gluon distribution used as input in (§) and @. In general, 
it is assumed that the gluon distribution is described by a parametrization 
of the parton distributions (for example: GRV, MRS, CTEQ). In this case 
the unitarity corrections in the gluon distribution are not included explicitly. 
However, calculating the unitarity corrections using the AGL approach we 
obtain that they imply large corrections in the behavior of the gluon distribu- 
tion in the HERA kinematical region, and therefore cannot be disregarded. 
Therefore, we should consider the solution from Eq. (|1|) as input in the Eqs. 
(§) and (U) in order to accurately determine the behavior of F2 and its slope, 
i.e., we should consider the unitarity corrections in the quark and the gluon 
sectors (quark+ gluon sector). 

The expressions (|J) , @ and (U) are correct in the double leading logarith- 
mic approximation (DLLA). As shown in |TI[] the DLLA does not work quite 
well in the whole accessible kinematic region (Q 2 > OAGeV 2 and x > 10~ 6 ). 
Consequently, a more realistic approach must be considered to calculate the 



observables in the HERA kinematical region. In [11] the subtraction of the 



Born term and the addition of the GRV parametrization |T7|] were proposed 
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to the F 2 case. In this case we have 



F 2 (x,Q 2 ) =F 2 (x,Q 2 )[Eq. (3)]-F 2 (x,Q 2 )[Born] + F 2 (x,Q 2 )[GRV] , (5) 
where the Born term is the first term in the expansion in K q of the equations 
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18| for more details). Here we apply this procedure for the F 2 slope. 



In this case 

dJ ^f ~ <41 " ^ + *g£Wl , (a) 

where the Born term is the first term in the expansion in K q of the equation 
The last term is associated with the traditional DGLAP framework, 
which at small values of x predicts 

dlogQ* ~ 9vr Jo azW) \-z 9 \\-z M J ' [ ) 

where a s (Q 2 ) is the running coupling constant and the splitting function 
P qg {x) gives the probability to find a quark with momentum fraction x inside 
a gluon. This equation describes the scaling violations of the proton structure 
function in terms of the gluon distribution. We use the GRV parametrization 
as input in the expression (0). 

In Fig. |2] we compare our predictions for unitarity corrections in the F 2 
structure function and the HI data |19| as a function of In(-) at different 



virtualities and some values of the radius. We see clearly that the unitarity 
corrections strongly increase at small values of the radius R. Our goal is not 
a best fit, but to eliminate some values of radius comparing the predictions 
of the AGL approach and HERA data. The choice R 2 = 1.5 GeV 2 does not 
describe the data, i.e. the data discard the possibility of very large SC in the 
HERA kinematic region. However, there are still a large range (5 < R 2 < 12) 
of possible values for the radius which reasonably describe the F 2 data. To 
discriminate between these possibilities we must consider the behavior of the 
F 2 slope. 

In Fig. [3] (a) we present our results for dF ^Q2 ^ considering initially the 
unitarity corrections only in the quark sector, i.e. using the GRV param- 
eterization for the gluon distribution as input in Eq. (f|). The ZEUS data 
points [2(J correspond to different x and Q 2 value. The (x, Q 2 ) points are 



averaged values obtained from each of the experimental data distribution 
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bins. Only the data points with < Q 2 > > 0.52 GeV 2 and x < 10 _1 were 
used here. Our results show that the fit to the data occurs at small values 
of R 2 , which are discarded by the F 2 data. Therefore, in agreement with our 
previous conclusions, we must consider the unitarity corrections to the gluon 
distribution to describe consistently the F 2 and dF ^ g q2 - data. In Fig. [3| (b) 

we present our results for dF ^ g q2 ~ considering the SC in the gluon and quark 
sectors for different values of R 2 , calculated using the AGL approach. The 
best result occurs for R 2 = 5GeV~ 2 , which also describes the F 2 data. 

The value for the squared radius R 2 = 5GeV~ 2 obtained in our analy- 
sis agrees with the estimates obtained using the HERA data on diffractive 
photoproduction of J/^f meson [21], 2^]. Indeed, the experimental values 



for the slope are B d = 4GeV~ 2 and B in = 1.66 GeV~ 2 and the cross sec- 
tion for J/ty diffractive production with and without photon dissociation are 
equal. Neglecting the t dependence of the pomeron-vector meson coupling 
the value of R 2 can be estimated |I5|. It turns out that R 2 ~ 5GeV~ 2 , i.e., 



approximately 2 times smaller than the radius of the proton. 

Some additional comments are in order. The unitarity corrections to F 2 
and its slope may also be analysed using a two radii model for the proton 
[ f23|j . Such analysis is motivated by the large difference between the measured 
slopes in elastic and inelastic diffractive leptoproduction of vector mesons in 
DIS. An analysis using the two radii model for the proton is not a goal of 
this paper, since a final conclusion on the correct model is still under debate. 

The AGL approach describe the ZEUS data, as well as the DGLAP evo- 
lution equations using modified parton distributions. Recently, the MRST 



24] and GRV group |25[] have proposed a different set of parton parametriza- 
tions which considers an initial 'valence-like' gluon distribution. In Fig. ^ 
we present the predictions of the DGLAP dynamics using the GRV(98) pa- 
rameterization as input. We can see that this parametrization allows to 
describe the F 2 slope data without an unconventional effect. This occurs 
because there is a large freedom in the initial parton distributions and the 
initial virtuality used in these parametrization, demonstrated by the large 
difference between the predictions obtained using the GRV(94) or GRV(98) 
parametrization. In our analysis we assume that the unitarity corrections 
are present in the HERA kinematical region, mainly in the F 2 slope, which 
is directly dependent of the behavior of the gluon distribution. Therefore, 
we consider that the fail of the DGLAP evolution equation plus GRV(94) to 
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describe the ZEUS data as an evidence of the high density effects and the 
possibility of description of the data using new parametrizations as a way to 
hidden these effects. In Fig. ^ we also show that if the GRV(98) parameter- 
ization is used as input in the calculations of the high density effects in the 
HERA kinematical region, the ZEUS data cannot be described. 

A comment related to the F 2 slope HERA data is important. The ZEUS 



data |2(J were obtained in a limited region of the phase space. Basically, 
in these data a value for the F 2 slope is given for a pair of values of x and 
Q 2 , i.e. if we plot the F 2 slope for fixed Q 2 as a function of x we have only 
one data point in the graphic. Recently, the HI collaboration has presented 



a preliminary set of data for the F 2 slope [26| obtained in a large region 
of the phase space. The main point is that these new data allows us an 
analysis of the behavior of the F 2 slope as a function of x at fixed Q 2 . In 
Fig. we present the comparison between the predictions of the DGLAP 
approach, using the GRV(94) or GRV(98) parametrization as input, and the 
AGL approach with the preliminary HI data [extracted from the Fig. 13 of 
[p6|] 1. We can see that, similarly to observed in the ZEUS data, the DGLAP 
+ GRV(94) prediction cannot describe the data, while the AGL approach 
describe very well this set of data. We believe that our conclusions are not 
modified if these new data in a large phase space are included in the analysis 
made in this letter. 

In this paper we have assumed that the unitarity corrections (high density 
effects) are present in the HERA kinematical region and believe that only a 
comprehensive analysis of distinct observables (Fl, F 2 , dF ^Q2 ^ ) will allow a 
more careful evaluation of the unitarity corrections at small x [|i~8j 27\ . The 
main conclusion of this paper is that the analysis of the F 2 and dF ^f q2 
data using the AGL approach implies that the gluons are not distributed 
uniformly in the whole proton disc, but behave as concentrated in smaller 
regions. This conclusion motivates an analysis of the jet production, which 
probes smaller regions within the proton, in a calculation that includes the 
high density effects. 
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Figure Captions 



Fig. |} First order contribution to the unitarity corrections. In (a) these 
corrections are controlled by the proton radius, while in (b) by the constituent 
radius. 



Fig. |[ The F2 structure function as a function of the variable In(-) for 
different virtualities and radii. Only the unitarity corrections in the quark 
sector are considered. Data from HI |19 . 



Fig. [| The F2 slope as a function of the variable x for different radii, (a) 
Only the unitarity corrections in the quark sector are considered, (b) The 
unitarity corrections in the gluon-quark sector are considered. Data from 
ZEUS [EOl. The data points correspond to a different x and Q 2 value. 



Fig. |j: Comparison between the DGLAP and Glauber-Mueller (GM) 
predictions for the behavior of the F2 slope using as input in the calculations 
the GRV(94) or GRV(98) parameterizations. Data from ZEUS [[2(J. The 
data points correspond to a different x and Q 2 value. 



Fig. ^ Comparison between the Glauber-Mueller (GM) prediction and 
DGLAP using as input in the calculations the GRV(94) or GRV(98) param- 
eterizations, for the behavior of the F2 slope. Preliminary data from HI 
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